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Very recently we have developed a dynamic network model for eco-systems that achieved "uni- 
fication" of "micro" and "macro"-evolution. We now propose an extension of our model so as to 
stabilize the eco-system and describe speciation in a more realistic manner. 
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I. INTRODUCTION 



It is now widely believed that functional networks abound in biological systems from molecular level (e.g., 

genetic and metabolic networks) to the levels of cells, organisms and species. An eco-system may be viewed as 
a functional network of species where the (directed) links of the network represent the inter-species interactions 
0. These interactions include not only prey-predator interactions but also competitions between species, for the 
same resources, as well as possible cooperations. Each species is also a network of individual organisms; the intra- 
species interactions can be cooperative or competitive. The problems of evolutionary ecology have been investigated 
extensively from theperspective of statistical physics |1, S lo] ■ 

In recent papers 0, 0| we have developed a model from the perspective of dynamic networks, by incorporating 
both inter-species and intra-species interactions, for studying some generic features of the biological evolution of eco- 
systems. In this paper, following a brief review of this model and its main results, we propose an extension of the 
model so as to capture speciation in a more realistic manner. 



II. EARLIER MODELS AND THEIR LIMITATIONS 



For convenience, most of the earlier theories focussed attention either on "micro"- evolution on ecological time 
scales (e.g., annual variations in the populations of species, mortality rates, etc.) or on "macro" -evolution that is 
most prominent on geological times scales (e.g., the phenomena of speciation or extinction). The ecological models 
that describe population dynamics using, for example, Lotka-Volterra type equations usually ignore macro- 
evolutionary changes in the eco-system. On the other hand, most of the macro-evolutionary models 0, 0, 0| 
do not explicitly explore the ageing and age-distributions of the populations of various species in the system. The 
models of ageing and dynamics of age-distributed populations ,14] , usually, focus on only one single species and do 
not incorporate the inter-species interactions that are, however, crucially important for their extinctions. 

But, in reality, evolution is a continuous process. When monitored at short intervals and over a not-too-long period 
of time, the ecological processes dominate the visible changes in the eco-system. However, if the same system is 
watched over sufficiently long period of time, the macro-evolutionary changes cannot be missed. Therefore, it is 
desirable to have a "unified" theoretical model that would unfold the natural continuous process while simulated on 
a computer. Very recently, attempts have been made b y se veral group s to merge population dynamics and "macro" - 
evolution within a single mathematical framework 15 Il5l IT^ HtL Hq . These efforts may have been made feasible, at 
least in part, because of the availability of fast computers. In our models 0,13 we have achieved not only the merger 
of population dynamics and "macro" -evolutionary processes but also detailed description of birth, ageing and death 
of individual organism so that the age-distributions in the populations of different species can also be monitored. 

Moreover, for the study of population dynamics of t he sp ecies in the eco-system one needs a model of the food web, 
a graphic description of prey- predator relations [j, ll^, ■ More precisely, a food web is a directed graph where 
each node is labelled by a species' name and each directed link indicates the direction of flow of nutrient (i.e., from 
a prey to one of its predators). In contrast to most of the contemporary models published in the physics literature 
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(ref.|23| is one of the few exceptions), we incorporate the trophic level structures of real food webs through a generic 
hierarchical model (81 . Besides, in order to capture ecological and evolutionary processes within the same theoretical 
framework, we allow the food web to evolve slowly with time |2^ . 

III. THE NETWORK MODEL: COMPONENT AND MOTIVATIONS 
A. Architecture of the network 

We model the eco-system as a dynamic network each node of which represents a niche that can be occupied by 
at most one species at a time. The network considered in our earliest formulation of the " unified" model jTj can be 
schematically represented by a random network. In a subsequent paper Q we replaced the random architecture by a 
generic hierarchical one, where niches are arranged in different trophic levels, with biologically realistic inter-species 
interactions. The hierarchical architecture helps us in capturing a well known fact that in the normal ecosystems the 
higher is the trophic level the fewer are the number of species. 

We assume only one single species at the highest level £ — 1. There are m^~^ levels at the ^-th level where m is 
a positive integer. The allowed range of £ is 1 < < imax- At any arbitrary instant of time t the model consists of 
N(t) species each of which occupies one of the nodes of the dynamic network. The total number of species cannot 
exceed Nmax = (m^™"^ — l)/{m — 1), the total number of nodes. Our model allows N[t) to fluctuate with time over 
the range I < N{t) < N,nax- The population (i.e., the total number of organisms) of a given species, say, i, at any 
arbitrary instant of time t is given by ni(t) < Umax- Thus, the total number of organisms n{t) at time t is given 

by n{t) = Yld-^t' Note that Imaxi (and, therefore, N^nax) and Umax are time-independent parameters in the 

model. 

The network itself evolves slowly over sufficiently long time scales. For example, random genetic mutations are 
captured by implementing random tinkering of some of the intra-node characteristics which will be introduced in 
the next subsection. The inter-node interactions change slowly to capture adaptive co-evolution of species thereby 
altering the graph that represents the network. Even the occupants of the nodes can change with time because, 
as the eco-system evolves, the populations of some species would drop to zero, indicating their extinction, and the 
corresponding nodes would be slowly re-occupied by new species through the process of speciation. 

B. Intra-node characteristics, mutations and dynamics 

The faster dynamics within each node captures "micro" -evolution, i.e., the birth, growth (ageing) and natural 
death of the individual organisms. For simplicity, we assume the reproductions to be asexual. An arbitrary species i 
is collectively characterized by 0: 

(i) the ijiinimum reproduction age ^|.ep(z), 

(ii) the birth rate M{i), 

(iii) the maximum possible age Xmaxii) — 100x2^^^^^/^ that depends only on the traophic level occupied by the species. 

An individual of the species i can reproduce only after attaining the age Xrep{i)- Whenever an organism of this 
species gives birth to offsprings, M{i) of these are born simultaneously. None of the individuals of this species can 
live longer than Xmax{i), even if an individual manages to escape its predators. The explicit form of X^ax assumed 
above is intended to mimic the fact that the species at the higher trophic level usually have higher lifespan. 

During each time step, because of random genetic mutations, Xrep and M independently increase or decrease by 
unity, with equal probability, Pmut ■ Xrep is not allowed to exceed a predetermined (large) positive integer while M is 
restricted to remain positive. 

The mtra-species competitions among the organisms of the same species for limited availability of resources, other 
than food, imposes an upper limit n^nax of the allowed population of each species. In order to capture this requirement, 
we assume the time- dependent probability pf,{i,a) (of individual a in species i) of giving birth per unit time to be 
a product of two factors; one of the factors is a standard Verhulst factor 1 — ni/n„iax whereas the other factor 
takes into account the age-dependence of Pb{i^ oi). For an organisms, this second factor becomes non-zero only when 
it attains the corresponding minimum reproductive age X^ep] this non-zero factor is assumed to be of the form 
{Xmax — X) I (Xmax — A^rep)- Notc that the latter factor is unity sX X = X^ep, and thereafter decreases, linearly 
with age, to zero at Xmax- Thus in the limit of vanishingly small population, i.e., Ui — > 0, we have pb{i,a) ^ 1 if 
X{i,a) = Xrep{i) and, thereafter, pi, decreases linearly |24| as the organism grows older. However, pi,{i,a;t) — > as 
ni{t) njnax, irrespective of the age of the individual organism a |25|. Occasionally, random mutations mentioned 
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above can lead to anomalous situations where some organisms may have Xrep > Xmax', for such organisms Pb{i, a; t) = 
for all t and these fail to reproduce during their entire life time. 

Similarly, we assume the probability pd of "natural" death (due to ageing) to be a constant pd — exTp[—r{Xmax — 
Xrep)/M], where r is a small fraction, so long as X < Xrep- However, for X > Xrep, the probability of natural death 
increases exponentially (following Gompertz law) as pd — exp[—r{Xmax — X)/M]. Note that, for a given Xmax and 
Xrep, the larger is the M the higher is the pd for any age X. Therefore, in order maximize reproductive success, each 
species has a tendency to increase M for giving birth to larger number of offsprings whereas the higher mortality for 
higher M opposes this tendency 

Because of the natural death mentioned above and, more importantly, prey-predator interactions (to be described 
in the next subsection), the populations of some species may fall to zero. The nodes left empty by such extinct species 
are then re-filled by new species. In order to capture this process of speciation, all the empty nodes in a trophic level 
of the network are re-filled by random mutants of one common ancestor which is picked up randomly from among 
the non-extinct species at the same trophic level. The subsequent accumulation of random mutations over sufficiently 
long time leads to the divergence of the genomes of the parent and daughter species that is an essential feature of 
speciation. 

However, occasionally, all the niches at a level may lie vacant. Under such circumstances, all these vacant nodes are 
to be filled by a mutant of the non-extinct species occupying the closest lower level. In our computer simulations, the 
search for this non-extinct species is carried out in steps, if even the lower level is also completely empty, the search 
for survivor shifts to the next lower level and the process continues till the lowest level is reached. If all the nodes, 
starting from the lowest, upto a certain level ever fall vacant, then, no new speciation takes place and the starvation 
deaths of the species propagate up the layers ending, finally, with the collapse of the entire eco-system. 

C. Inter-node interactions and dynamics 

The interaction between any two species i, k that occupy two adjacent trophic levels is given by Jik- The sign of J^^ 
gives the direction of trophic flow, i.e. it is -1-1 Hi eats k and it is —1 if fc eats i. In the absence of any prey-predator 
interaction between the species i and k, Jik = 0. For simplicity, we assume the absolute value (magnitudes) of all the 
non-vanishing interactions to be unity. Note that although there is no direct interaction between species at the same 
trophic level in our model, they can compete, albeit indirectly, with each other for the same food resources available 
in the form of prey at the next lower trophic level. 

The J account not only for the mfer-species interactions but also mfra-species competitions for food. Let be 
the number of all prey individuals for species i on the lower trophic level, and S~' be m times the number of all 
predator individuals on the higher trophic level. Because of the larger body size of the predators |23| , we assume that 
a predator eats m prey per time interval. Then, gives the available food for species i, and S~ is the contribution 
of species i to the available food for all predators on the next higher level. If the available food is less than the 
requirement, then some organisms of the species i will die of starvation, even if none of them is killed by any predator. 

If Ui — is larger than S~ then food shortage will be the dominant cause of premature death of a fraction of the 
existing population of the species i. On the other hand, if > Ui — , then a fraction of the existing population 
will be wiped out primarily by the predators. 

It is well known that each species tries to minimize predators but, at the same time, looks for new food resources. 
In order to capture this, at each time step, each of the species in our model, with the probability Pmut, re-adjusts a 
link J from one of its predators and another to one of its potential preys . If the link Jij to the species i from 
a higher level species is non-zero, it is reassigned a new value = Jji = 0. On the other hand, if the link Jik to a 
species i from a lower level species k was zero, the new values assigned are Jik = 1, Jki = ^1- 

D. Initial conditions and update rules 

The requirements of computational resources increase exponentially with increasing £rnax- Therefore, in almost all 
our simulations we chose £max = 5, although we verified that the qualitative features of the data were similar in case 
of smaller imax- In our simulations, we always began with a random initial condition where M = 1 for all species. 
Since larger species occupy the higher tropic levels and are expected to live longer than those at lower levels, we 
assinged Xmax = 100,71,50,35,25 to the species at level £ = 1,2,3,4,5, respectively. Initially, Xrep was assigned 
randomly between 1 and Xmax, the population randomly between 1 and nmax/'2- The ages of the individuals in the 
initial state varied randomly between 1 and the Xmax of the corresponding species. 

The state of the system is updated in discrete time steps as follows: 
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64000 to 64 samples, t=10M,10'^5,10'^6,10'^7 (n_max=100) and 10'^5 (n_max=1000, 320 samples) 
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FIG. 1: Log- log plots of the distributions of the lifetimes of the species. The common parameters for all the curves are 
m = 2,£ = 5 (i.e. Nmax = 31), psp = 0.1,Pmut = 0.001, C = 0.05, r = 0.05. The symbols +, x and * correspond to 
t = 10'',10^10'',10^ averaged over 64000, 6200, 640 and 64 systems respectively. The upper curves axe all for nmax = 100 
whereas the lower curve is for rimax = 1000. The line with slope —2 corresponds to a power law distribution that has been 
predicted by many theories. 

Step I- Birth: Each individual organism a {a = 1, ni{t)) of the species i {i = 1, 2, ...N{t)) is allowed to give birth to 
M{i; i) offsprings at every time step t with probability (per unit time) Pb{i^ a; t) the explicit form of which has been 
mentioned above. 

Step II- Natural death: At any arbitrary time step t the probability (per unit time) of "natural" death (due to ageing) 
of an individual organism a of species i is Pd{i, a; t). 

Step III- Mutation: With probability Pmut per unit time, mutations of intra-node characteristics and the interactions 
J are implemented. 

Step IV- Starvation death and killing by prey: 
At every time step t, in addition to the natural death due to ageing, a further reduction of the population by 

C max(5r, m-Sf) (1) 

is implemented where ni{t) is the population of the species i that survives after the natural death step above. C is a 
constant of proportionality. If implementation of these steps makes Ui < 0, species i becomes extinct. 
Step V- Speciation: The nodes left empty by extinction are re-filled by new species, with probability Psp following the 
algorithm for speciation mentioned above. 

The longest runs in our computer simulations were continued upto 10^ time steps. If each time step in our model 
is assumed to correspond to a real time of the order of one year, then the total time for which we have monitored our 
model eco-system, is comparable to real geological times scales. 
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Force of mortality; same runs as for Fig.1 ; Gompertz law corresponds to straight line 
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FIG. 2: Semi-log plot of the mortality function against age of the individual organisms at the highest trophic level. The same 
symbols in fig^and fig |5| correspond to the same set of parameter values. 

IV. RESULTS 
A. Lifetime distributions 

The average distributions of the hfetimes of the species are plotted in figQ] It is not possible to fit a straight line 
through the data over the entire range of hfetimes ahhough only a limited regime is consistent with a power-law with 
the effective exponent 2, which has been predicted by several models of "macro" -evolution 0,0. The overlap of the 
curves for different simulation times establishes that our simulations have reached the asymptotic regime where effects 
of initial conditions have been completely washed out. Because of the various known limitations of the available fossil 
data, it is questionable whether real extinctions follow power laws and, if so, over how many orders of magnitude. 

B. Distributions of Species Characteristics 

We define the mortality rate by the relation —d \n{survivors) / d{age). The mortality rate extracted from the raw 
data is plotted against age in fig[21 The shape of this curve is consistent with the usual census data that show a 
minimum in childhood and exponential increase in adults. 

During the early stages of the macro-evolution, the distributions of M broadens. But, with further passage of 
time it shrinks and reaches a stationary form where the largest M is around 30 (for the parameter set used in our 
simulations) . This is consistent with the view that organisms have a choice of either faster reproduction and shorter 
life or slower reproduction and longer life span. 
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Same runs as for lifetimes; squares give < 1 000 indiv., otiier symbols < 1 00 indiv, + for t = 1 0'^4 only 
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FIG. 3: Semi-log plot of the distribution of Al. The same symbols in figQand fig|3 correspond to the same set of parameter 
values. 



C. Collapse of fragile ecosystems 



One crucial effect of the generalization of the rule for speciation proposed in this paper is that the eco-systems are 
now much more stabilized than in our earlier papers 0j 13 ■ In our earlier work, reported in ref. we allowed an empty 
node to be re-occupied by one of the non-extinct species from only the same level where the extinct node was located. 
If, by chance, all the nodes of one level fell vacant at some stage, no speciation could take place; this would trigger an 
avalanche of extinctions that would propagate upward in the food chain because of starvation of the organisms and 
eventually lead to a collapse of the entire eco-system. In contrast, in the current version of our model, speciation is 
allowed to take place from lower levels if all the nodes in an entire level become vacant. This reduces the possibility 
of collapse of the eco-system but does not rule it out completely. Work on further stabilization of the eco-system with 
variable £max is in progress and will be reported elsewhere in a future publication (29i |. 



V. SUMMARY AND CONCLUSION 



In summary, we have extended our dynamic network model for unified description of micro- and macro-evolution. 
A majority of the main characteristics of the model are emergent properties of the self-organizing dynamics of the 
system. The main effect of the extension of the speciation dynamics proposed in this paper is that even after large 
avalanches of extinctions the eco-system can recover by speciation and bio-diversification starting from the surviving 
species at the low trophic levels. Consequently, the complete collapse of the eco-system becomes extremely rare and 
the distribution of the litter size M became stationary much faster. 
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